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Abstract 

It is shown that the critical threshold for percolation of the overlapping 
strings exchanged in heavy ion collisions can naturally explain the sharp 
strong suppresion of J / ip shown by the experimental data on central Pb-Pb 
collisions, which does not occur in central O-U and S-U collisions. 
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The NA50 collaboration (0) has reported a strong suppression of J ftp production 
in Pb-Pb central collisions at 158 AGeV/c per nucleon. The suppression is much 
stronger than the expected one due to J / ip absorption corresponding to a cross section 
of 6.3 mb, by which the NA38 data for O-U and S-U central collisions ([£3, H|) and the 
hadron-nucleus data can be explained. The NA50 data show a clear deviation from the 
previous situation (0). The J/ip suppression in peripheral Pb-Pb collisions is similar 
to the one corresponding to central S-U collisions, but a sharp enhancement occurs as 
the centrality of the Pb-Pb collisions increases. 

In this paper we draw attention to the fact that the continuum percolation of colour 
strings can naturally describe the sharp difference in the J/ip suppression at present 
energies between O-U, S-U and peripheral Pb-Pb collisions on the one side and central 
Pb-Pb collisions on the other side. Predictions for RHIC and LHC energies are given. 

The continuum percolation of colour strings takes place when the density of strings 
rises above a threshold, which can be calculated on geometrical grounds. In this picture, 
the region where several strings fuse can be considered a droplet of a non-thermalized 
Quark Gluon Plasma, in which the J /if} is suppressed as predicted by Matsui and Satz 
(||). Percolation means that these droplets overlap and the Quark Gluon Plasma 
domain becomes comparable to the nuclear size. 

In many models of hadronic collisions (|]6]|-||11||), colour strings are exchanged be- 
tween projectile and target. The number of strings grows with the energy and with 
the number of nucleons of the participant nuclei. When the density of strings becomes 
high the string colour fields begin to overlap and eventually individual strings may fuse 
(|12|-|]16|), forming a new string which has a higher colour charge at its ends, corre- 
sponding to the summation of the colour charges located at the ends of the original 
strings. The new strings break into hadrons according to their higher colour. As a 
result, heavy flavour is produced more efficiently and there is a reduction of the total 
multiplicity (|13|]). Also, as the energy-momenta of the original strings are summed to 
obtain the energy-momentum of the resulting string, the fragmentation of the latter 
can produce some particles outside the kinematical limits of nucleon-nucleon collisions 
if the original strings come from different nucleons ( [Jl7|, |18| ) . The fusion of strings has 
been incorporated in several Monte Carlo codes. In particular, in the Quark Gluon 
String Model (QGSM) it is assumed that strings fuse when their transverse positions 
come within a certain interaction area a ([0). The value of a is determined to repro- 
duce A rapidity distributions in S-S and S-Ag central collisions at pi a b= 200 GeV/c 
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and in Pb-Pb central collisions at pi a b= 158 GeV/c. 

Cascade reactions like 7r + p— >K°A, 7r°n— ►i^A, 7r + A-^K + p andpA— >tt + K° also con- 
tribute to the A rapidity distribution but their effects are smaller than the ones due to 
string fusion, generating uncertainties in the value of a of around 10%. From the value 
of a, the radius r of the transverse dimension of the string can be obtained, a = 2iir 2 
(jl9|). In our code only fusion of two strings is considered, so the obtained r-value, 
r=0.36 fm, is an effective one, somewhat larger than the real transverse radius of the 
string. Denoting by Nj the number of strings which fuse into j-fold strings and N 2 
and r e jf the number of all fused strings and the effective transverse size of the string, 
respectively, we will have 

2N' 2 ixrl ff = Y.N^r\ (1) 

3=2 

N' 2 = J2N r (2) 

J'=2 

The upper limit of the sum in (1) is determined by the constraint (2). The values 
of N 2 and r 2 tf were fixed in our calculation by comparing the results of the string 
fusion model with the experimental data on A production in central S-S collisions at 
y^s =19.4 AGeV. Computing Nj in our Monte Carlo code we obtain from (1) the value 
r=0.2 fm both for Pb-Pb and S-Ag collisions. 

In nucleus-nucleus collisions many strings are exchanged. In impact parameter 
space these strings are seen as circles inside the total collision area. As the number 
of strings increases, more strings overlap. Several fused strings can be considered as 
a domain of a non-thermalized Quark Gluon Plasma. Following the arguments of 
Matsui and Satz (||) the J/ip can not be formed inside this domain. Also, the J/ip 
will be destroyed by interaction with these fused strings. Above a critical density of 
strings percolation occurs, so that paths of overlapping circles are formed through the 
whole collision area. Along these paths the medium behaves like a colour conductor. 
The percolation gives rise to the formation of Quark Gluon Plasma on a nuclear scale. 
The phenomenon of continuum percolation is well known (|[20||). It explains hopping 
conduction in doped semiconductors and other important physical processes ( |[21]| ). The 
percolation threshold r] c is related to the critical density of circles n c by the expression 

Vc = vrr 2 n c . (3) 
T) c has been computed using Monte Carlo simulation, direct-connectedness expansion 
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and other different methods. All the results are in the range r] c = 1.12 — 1.175 (||22j|- 
26|). Taking the above mentioned value of r, these values imply 



n c = 8.9 — 9.3 strings/fm 2 . 



(4) 



One may introduce a hard core to model a repulsive interaction between the circles, 
or to substitute circles by squares. The percolation threshold r\ c is only slightly reduced 
in these cases. This enhances the confidence in its value and the application to our 
case where we do not know the dynamics of the interaction among strings. 

In Table 1 the number of strings exchanged for central p-p, S-S, S-U and Pb-Pb 
collisions is shown together with their densities. It is seen that at SPS energies only 
the density reached in central Pb-Pb collisions is above the critical density. In Pb-Pb 
minimun bias collisions the average number of strings at SPS energies is 227, very 
similar to the value for central S-U collisions, so the density is lower than the critical 
one. 

The J I ip suppression experimentally observed follows the same pattern. The strong 
suppression is only observed in central Pb-Pb collisions. According to Table 1, a strong 
J/ip suppression is also expected in S-U collisions at RHIC energies and in S-S and 
S-U collisions at LHC energies. 

Recently (f27|) it has been assumed that the produced J/ij) is completely destroyed 
whenever the energy density exceeds a certain value and this energy density is taken 
proportional to the density of participants. The critical value is chosen to lie between 
the density of participants of central S-U collisions and Pb-Pb collisions. With this 
choice a good description of the experimental data is obtained. In our model the 
density of strings is proportional to the number of collisions, and we obtain similar 
quantitative results. However, in our approach the critical value is naturally explained 
on geometrical grounds. 

In Fig. 1 the distribution of strings fusing into sets of a given number of fused 
strings is shown for central S-U collisions at -\/i=19.4 AGeV and i/i=200 AGeV and 
also for central Pb-Pb collisions at a/s=19.4 AGeV. The first case is below and the 
second above the percolation threshold. It is seen that above the percolation threshold 
we can obtain many sets with a very high number of fused strings. 

Refering to tp' suppression the experimental data reveal the following features (U, 
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1) The ratio ip' /ip is constant in p-A collisions. 

2) ip' /ip decreases with centrality in S-U collisions. The decrease seems to stop at high 
centrality. 

3) ip' /ip is almost the same in central Pb-Pb and S-U collisions. 

The first two features of experimental data can be explained by absorption and 
interaction with comovers (|[31||-p3|). Taking equal absorption cross section a(J/ip) = 
a{ip') ~ 4.2 mb, the hadron-nucleus behaviour can be explained since no interaction 
with the produced particles (f33j) is assumed. In nucleus-nucleus collisions low energy 
interactions of J / ip and ip with the hadrons produced in the collision break both the 
ip' and the J / ip but the ip cross section at low energy is much larger than that of J / ip 
(the threshold for breaking the ip is only 52 MeV and the one for J ftp is 640 MeV). 
This difference in the cross sections may be responsible for the different behaviour of 
ip' and J I ip suppression in central S-U collisions. In our picture this behaviour can be 
explained by noting that in S-U central collisions the average distance between strings 
is the order of 0.4 fm, larger than the size of J /ip (0.2 fm) but less than the size of the ip . 
Therefore one expects that ip interacts with the strings or with the particles produced 
by the strings with greater probability than J/ip. For central Pb-Pb collisions, with 
the density above the critical percolation threshold, no additional suppression of ip' 
relative to J/ip is expected, in agreement with the data. 

Also it is possible that the percolation process takes place among the produced 
resonances and particles instead of strings (fl3"4"|). The two cases can be distinguished 
by studying the behaviour of long range correlations and measuring forward-backward 
correlations (p5*|). 

The percolation of strings can be considered as a smooth way to Quark Gluon 
Plasma. Around percolation threshold, strong fluctuations in the number of strings 
with a given colour should appear. This will produce large fluctuations in a number 
of different observables, like strangeness, in an event by event analysis. Also a large 
number of Or (confirmed by the experimental data (|36[])) and a copious production 
of hadronic particles with \xp\ much larger than 1, outside the kinematical nucleon- 
nucleon limits, may serve as clear signatures. The latter would also distinguish our 
picture from the percolation of resonances and particles. 
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Table captions 

Table 1. Number of strings (upper numbers) and their densities (fm~ 2 ) (lower num- 
bers) in central p-p, S-S, S-U and Pb-Pb collisions at SPS, RHIC and LHC energies. 
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Figure captions 

Figure 1. Percentage of the total number of strings exchanged in the collision which 
goes into sets of a given number of fused strings, for central S-U collisions at ^=19.4 
AGeV (dashed line) and ^=200 AGeV (solid line) and for central Pb-Pb collisions 
at y/s=19A AGeV (dotted line). The number 10 in the horizontal axis indicates sets 
of 10 or more strings. 
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Table 1 
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